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Abstract. We perform a first principles band calculation for a group of quasi- two- 
dimensional organic conductors ^-(BDA-TTP)2MF6 (M=P, As, Sb, Ta). The ab-initio 
calculation shows that the density of states (DOS) is correlated with the band width 
of singly occupied (highest) molecular orbital (SOMO), while it is not necessarily 
correlated with the unit cell volume. The direction of the major axis of the cross 
section of the Fermi surface lies in the F-B direction, which differs from that obtained 
by the extended Hiickel calculation. Then, we construct a tight-binding model which 
accurately reproduces the ab-initio band structure. The obtained transfer energies give 
smaller dimerization than in the extended Hiickel band. As for the difference of the 
anisotropy of the Fermi surface, the transfer energies along the inter-stacking direction 
are smaller than those obtained in the extended Hiickel calculation. Assuming spin- 
fluctuation-mediated superconductivity, we apply random phase approximation (RPA) 
to a two-band Hubbard model. This two-band Hubbard model is composed of the 
tight-binding model derived from the first principles band structure and an on-site 
(intra- molecule) repulsive interaction taken as a variable parameter. The obtained 
superconducting gap changes sign four times along the Fermi surface like in a d- 
wave gap, and the nodal direction is different from that obtained in the extended 
Hiickel model. Anion dependence of Tc is qualitatively consistent with the experimental 
observation. 
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1. Introduction 

Ever since tlieir discovery, pliysics and cliemistry of organic superconductors have 
attracted much attention [1]. One of the interesting features of these superconductors 
is the possibihty of unconventional pairing arising from electron correlation effects. 
From this point of view, there have been efforts to control the strength of the 
electron correlation via chemical modification of molecules, and this has lead to 
the discovery of several new superconductors [2J. Among those superconductors, 
here we focus on a group of organic conductors /3-(BDA-TTP)2MF6 (M =P, As, 
Sb, Ta), where BDA-TTP is an abbreviation of 2,5-bis(l,3-dithian-2-ylidene)-l,3,4,6- 
tetrathiapentalene. BDA-TTP molecule is derived from BDH-TTP by stereochemical 
modification, where BDH-TTP is an abbreviation of 2,5-bis(l,3-dithiolan-2-ylidene)- 
1,3,4,6-tetrathiapentalene. While BDH-TTP salts show stable metallic behavior 
suggesting weak electron correlation, the BDA-TTP salts exhibit stronger electron 
correlation, which is expected from the band width narrowing due to the chemical 
modification. In fact, /3- (BDA-TTP) 2MF6 exhibit superconductivity with Tc=7.5, 5.8, 
5.9K for M =Sb, As, P, respectively 0, while M =Ta does not superconduct jl]. The 
conducting layer is sandwiched by insulating layers of anion MFg as shown in figure 
[U^a), from which we may expect these salts to be quasi- two-dimensional. The alignment 
of the BDA-TTP donors in the conducting layer takes the /3-type configuration shown 
in figure [D^b). Measurements of Shubnikov-de Haas (SdH) oscillation and angular- 






Figure 1. Lattice structure of /3-(BDA-TTP)2SbF6: (a) the side view and (b) the 
top view of the conducting BDA-TTP layer. 



dependent magneto-resistance oscillation (AMRO) for /3-(BDA-TTP)2SbF6 have shown 
a two-dimensional Fermi surface [5] , which is consistent with the tight-binding band 
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structure [5] obtained from the extended Hiickel calculation [B]. However, a very recent 
AMRO measurement shows that the direction of the anisotropy of the Fermi surface 
differs from the previous AMRO measurement, and suggests the direction of the major 
axis of the cross section lies in the F-B direction of the Brillouin zone (see figure [2]) 
[7]. It is worth mentioning that the effective cyclotron mass takes a large value of 
ml /mo = 12. 4 ±1.1, which indicates that these materials are strongly correlated electron 
systems, although the normal state exhibits a metallic nature. 

Measurements of the upper critical field {Hc2) of /3-(BDA-TTP)2SbF6 have shown 
spin-singlet pairing [SI [9] . The temperature dependence of the electronic specific heat 
has indicated presence of a nodal structure in the superconducting order parameter 
[8], and a scanning tunneling microscope (STM) experiment also suggests presence of 
anisotropy in the gap, where the nodes (or the gap minima) lie in the ka ± kc direction 
|10] . On the other hand, a recent study of the in-plane anisotropy of Hc2 shows that 
the Hc2 maxima appears when the magnetic field is applied in the ka ± kc direction, i.e., 
the nodal direction indicated in the STM measurement [H]. This is puzzling because 
the Hc2 maxima usually appears when the field is applied in the antinodal direction. A 
uniaxial pressure experiment for M =As and Sb shows that Tc once increases and takes 
a maximum under the compression parallel to the stacking direction [12] . 

Theoretically, several studies have adopted the models derived from extended 
Hiickel calculation, and analysis on superconductivity has also been performed assuming 
spin fluctuation mediated pairing [131 HI]- Nonoyama et al. [151 US] applied RPA to 
the two band model for M =As and Sb, while Suzuki et al. [17] applied the fluctuation 
exchange (FLEX) approximation to the original two-band model and the single-band 
dimer model for M =As and Sb, and discussed the pressure dependence. There, the 
(i-wave like gap function has been obtained, whose nodal direction is consistent with 
the STM measurement [10] , but not with the recent study of in-plane anisotropy of Hc2 

Thus, there are several controversies regarding the relation between experiments 
and theoretical studies based on the extended Hiickel model. In the present paper, we 
study theoretically the organic conductors /3-(BDA-TTP)2MF6 (M =P, As, Sb, Ta) 
based on a model constructed from ab-initio band calculation based on the density 
functional theory (DFT) exploiting the WIEN2K package ^S\- In fact, the importance 
of model construction of molecular solids based on ab-initio calculation has recently 
pointed out in several studies [191 EOl ET]. We point out the difference in the band 
structure and the Fermi surface between the ab-initio (non-empirical method) and 
the extended Hiickel calculations (semi-empirical method), and its consequence to the 
superconducting gap mediated by spin fluctuations. We also study the anion dependence 
of the band structure, and how it can affect the superconducting transition temperature. 
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2. Method 

2.1. ab-initio band calculation and model construction 

We perform an ab-initio band calculation using all-electron full potential linearized 
augmented plane- wave (LAPW) + local orbitals (lo) method within WIEN2K [18]. This 
implements the DFT with different possible approximation for the exchange correlation 
potentials. The exchange correlation potential is calculated using the generalized 
gradient approximation (GGA). 

To attain convergence in the eigenvalue calculation, the single-particle wave 
functions in the interstitial region are expanded by plane waves with a cut-off of 
i?MT-^max = 3.0, whcre -Rmt denotes the smallest muffin tin radius and -ft'max is the 
maximum value of K vector in the plane wave expansion. For /3-(BDA-TTP)2SbF6, 
for instance, the muffin-tin radii are taken as 1.74, 1.74, 1.62, 0.83, and 0.45 in atomic 
units (au) for Sb, F, S, C, and H, respectively. Thus -ft'max is 3/0.45~6.7, and the plane 
wave cutoff energy is 604.7 eV. Calculations are performed using 7x3x9 fc-points in the 
irreducible Brillouin zone. We adopt the lattice structure determined experimentally 
(including the coordinates of the H atoms) [3] , and we do not relax the atomic positions 
in the calculation. 

Having done the ab-initio calculation, we then construct a tight-binding model 
which accurately reproduces the ab-initio band structure. From figure [TJ^b), it can be 
seen that the conducting BDA-TTP layer has two BDA-TTP molecules in a unit cell, 
so we regard one molecule as a site, and consider a two-band (two sites per unit cell) 
tight-binding model to fit the ab-initio band structure. The tight-binding Hamiltonian 
we consider, Hq, is written in the form 

^0= E {tia-.j/BcLaCj^a + H-C.} , (l) 

where i and j are unit cell indices, a and (3 specifies the sites in a unit cell, c|„^ {ciaa ) 
is a creation (annihilation) operator with spin a at site a in the i-th unit cell, tia-.jp is 
the electron transfer energy between {i, a) site and (j, (3) site, and {ia : j(3) represents 
the summation over the bonds corresponding to the transfer. 
By Fourier transformation, equation ([T]) is rewritten as 

Ho= (^) (2) 

k,c7,a,l3 

where eai3 (k) is the site-indexed kinetic energy represented in fc-space. The band 
dispersion is given by diagonalizing Ea/s (k), 

£a/3 (fc) = E^"7 (^) ^^7 (^)^7 (^) ' (3) 

7 

where ^7 (A;) is the dispersion of the 7 band measured from the chemical potential, and 
da^y (k) is the unitary matrix that gives the transformation. 

In order to study unconventional superconductivity, we employ the two-band 
Hubbard model obtained by adding the on-site (intra-molecule) repulsive interaction 
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to the tight-binding model derived from the ab-initio band structure. The Hubbard 
Hamiltonian, H, is 

H = Ho + ^U nia^riiai (4) 

ia 

where U is the on-site interaction and riiaa is the number operator of the electron on 
a-site in the i-th unit cell. The on-site U in the present study is fixed at a certain 
value. Although this choice does not have a quantitative basis, taking other values 
of U does not affect qualitatively the conclusion of the present paper regarding the 
superconducting gap symmetry and the anion dependence of Tc. Note also that we 
ignore the off-site (inter-molecule) repulsive interactions. This is based on the result 
obtained in a previous theoretical work for the extended Hiickel model of the present 
materials [I5]. There it has been shown that the introduction of moderate inter-molecule 
interactions only slightly modifies the position of the nodes of the gap. Although the 
model adopted here is different, we believe that the effect of the inter-site interaction 
should be similar as far as the off-site interaction values are moderate. 



2.2. Random phase approximation and superconducting gap equation 

We apply RPA to the two-band Hubbard model given by equation (jl]) as follows. The 
bare susceptibility matrix in the site-representation is given by 

X%{q) = ^ dc^-r iP + Q) d}^ iP + Q) dpy (p) c?*y (p) 

P 7,7' 

where N is the total number of unit cells, and / (0 is the Fermi distribution function 
(this is the multi-site version of the Lindhard function). The indices a(3 on the left hand 
side means {a, /3)-element of the bare susceptibility matrix X^. Within RPA, the spin 
and charge susceptibilities are obtained as 

x'p = (^i -x^uy^ X^, (6) 

X=^= (/ + X°f/)"'x°, (7) 

respectively. X^^, X"^^, X° are the spin, charge, and bare susceptibility matrices, 
respectively, U is the on-site interaction matrix, and / is the unit matrix. In a two- 
band system, U, X'^p and X'^^ all become 2x2 matrices. The "spin susceptibility" 
(see figure O^a)) in the present paper is obtained as the larger eigenvalue of the 2x2 
spin susceptibility matrix. 

By using these susceptibility matrices, the pairing interaction v^^^s^'^^ for the spin- 
singlet state is given as 

ysinglet _ ^ ^ -UX'PU - -UX'^^U. (8) 
2 2 ^ ^ 

Using the pairing interaction given by equation (|8]), we solve the linearized gap equation 

to obtain the transition temperature Tc and the superconducting gap function for the 
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spin-singlet pairing state. The linearized gap equation within the weak-coupling theory 
is given by 

j,' a' 13' 77' 

X {k') rf*,^ (fc') d,3y {-k') d*p,y {-k') 

where A is the eigenvalue of the linearized gap equation, and ipai3 (k) is the (a, (3)- 
element of the gap function matrix. The transition temperature Tc is the temperature 
where A reaches unity. In a two-band system, ^/'^mgiet ^ become 2x2 matrices. 

Although RPA is quantitatively insufficient for evaluating the absolute value of T^, 
we expect this approach to be valid for (i) studying the form of the superconducting gap 
function assuming spin-fluctuation-mediated pairing, and (ii) qualitatively comparing 
the tendency toward superconductivity (or T^) between the different anions. We take 
128x128 fc-point meshes in the RPA calculation throughout the entire study. 



3. Results 



3.1. ab-initio band calculation 

We show in figure [2](a) the ab-initio band structure for /3-(BDA-TTP)2SbF6 at ambient 
pressure and room temperature. It can be seen that the singly occupied (highest) 
molecular orbital (SOMO) and the highest-occupied molecular orbital (HOMO) near the 
Fermi level (zero energy) are isolated from other bands, although there are two bands 
somewhat close to the HOMO band, which correspond to HOMO— 1 and HOMO— 2. 
Density of states (DOS) in figure [2](b) also shows the isolation of the HOMO and SOMO 
from the other bands. From these results, it can be considered that the SOMO and 
HOMO bands are the target bands to construct a low energy tight-binding model, 
figure [2](c) shows the Fermi surface of the ab-initio calculation, where the name of the 
fc-points are presented only on the ky^kb) = plane, e.g r=(0, 0, 0), Z=(0, 0, vr/c ) 
and B=(7r/a, 0, -vr/c ). We also show similar plots for /3-(BDA-TTP)2AsF6 at ambient 
pressure and room temperature in figure [31 The band structure and the Fermi surface 
are similar between the two salts. M =P and Ta also give similar results (not shown). 
The DOS also looks similar qualitatively, but its quantitative dependence on the anions 
will be discussed later. 

The Fermi surface is clearly cylindrical, reflecting the strong two dimensionality of 
this salt. The direction of the anisotropy of the Fermi surface differs from that found in 
the previous extended Hiickel calculation [31 HI [151 HZ] • The origin of this difference is 
intimately related to including/ excluding the 3d-orbital of the sulfur atoms of the TTP 
skeleton donors in the extended Hiickel calculation [22l [23l [24] . Actually, the shape of 
the Fermi surface given by the extended Hiickel calculation in which the 3d-orbital of 
the sulfur atoms is excluded [25] is similar to the present ab-initio result. Focusing on 
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Figure 2. (a) Calculated ab-initio band structure, (b) the density of states, and (c) 
the Fermi surface of /3-(BDA-TTP)2SbF6 at ambient pressure and room temperature, 
where X(Z)-axis corresponds to the a(c)-direction in the lattice structure. In figure 
(a), the red curve represent the ab-initio band dispersions and the blue solid curves 
gives the tight-binding fit. 

the comparison of the experiments with our results of the Fermi surface, the anisotropy 
of the Fermi surface differs from the previous experiment [5], but a very recent AMRO 
measurement shows that the long-axis of the Fermi surface ellipse is in the F-B direction, 
which is at least qualitatively consistent with the present result [7] . However, the size of 
the Fermi surface is not consistent between this experiment and the present calculation, 
and this remains as a puzzle. 

Now, one of the important issues we study in the present paper is the effect of 
substituting the anion. Here we consider the anion dependence of the band-width and 
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(a) . (b) 




Figure 3. Plots similar to Fi^for /3-(BDA-TTP)2AsF6 at ambient pressure and 
room temperature ; (a) Calculated ab-initio band structure, (b) the density of states, 
and (c) the Fermi surface. 

the DOS at the Fermi level for /3-(BDA-TTP)2MF6 (M=P, As, Sb, Ta) at ambient 
pressure and room temperature, and also for /3-(BDA-TTP)2SbF6 at 12K. As for the 
band width, we introduce H^sh and Wsomo, where the former is the energy difference 
between the top of SOMO and the bottom of HOMO, and the latter is the band width of 
SOMO alone, figure ID^a) shows that Wsn becomes smaller as the anion is enlarged from 
M =P to Sb, resulting in an increase of the unit-cell volume. The trend can naturally 
be understood as the effect of the chemical pressure. M =Ta, however, is an exception, 
where the increase of the unit-cell volume from M =Sb results in an increase in the 
band width. The relation between M =Sb and M =Sb(LT) shows that lowering the 
temperature results in an increases of W^sh due to thermal compression. We show in 
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Figure 4. (a) SOMO+HOMO band width, (b) SOMO band width, and (c) the 
density of states against the unit-cell volume. LT in parentheses for Sb means that 
the low temperature lattice parameters are used. The inset of (b) shows the schematic 
relation between SOMO and HOMO. 



figure m^b) Wsomo against the unit-cell volume. The correlation between PFsqmo and 
the volume is not monotonic in contrast to the volume dependence of Wsn, even within 
the pnictogens M =P, As and Sb. The reason for this lies in the small dimerization 
gap between the SOMO and HOMO bands. Even though PVsomo of M =As is larger 
than that of M =P, Wsr of As is smaller because the dimerization gap of As is smaller 
than that of P. A schematic figure is shown in the inset of figure IH^b), where the left 
(right) figure corresponds to M =P (As). We will come back to this point in the next 
subsection. The volume dependence of the DOS at the Fermi level, DOS(i?F), varies in 
accord with Hsomo, not H^sh, as seen in figure Ht^c). It is natural that DOS is inversely 
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related to Wsouo because this directly measures the width of the band that intersects 
the Fermi level, but the point here is that Wsomo is not directly correlated with W^n 
and thus the unit cell volume. 

3.2. Effective tight-binding model 

Figure |5] shows the effective tight-binding model adopted to fit the ab-initio band. In 
addition to the nearest-neighbor transfer energies (left part in figure [5l same notation 
is used as in previous studies), we need to introduce the next-nearest-neighbor transfer 
energies (right panel in figure |5]) to reproduce the ab-initio band more accurately. Note 
that the stacking direction of the BDA-TTP molecules is taken in the c -|- a-direction of 
the lattice structure in our study, as shown in figure [Dj^b). 




Figure 5. The tight-binding model adopted to fit the ab-initio band structure, where 
an ellipse corresponds to a BDA-TTP molecule. 

We show in figure |6]^a) the band structure of the two band tight-binding model for 
/3-(BDA-TTP)2SbF6 superposed to the ab-initio band structure. It can be seen that the 
tight-binding model nicely reproduces the ab-initio band structure, figure [6](b) shows 
the Fermi surface obtained from the model, which is essentially the same as the ab-initio 
result in figure [2t^c). 

Similar tight-binding fit is performed for all the anions, and the thereby determined 
transfer energies are summarized in Table [H The bottom two lines of Table [H represent 
the magnitude of the dimerization, tpi/tp2, in the stacking direction, and the relative 
transfer perpendicular to that direction (i.e. the strength of the two dimensionality), 
(tgi -\- tq2)/{tpi + tp2). We can now quantitatively discuss the effect of the dimerization 
on the unit-cell volume dependence of Wsouo and DOS(i?F)- Table [T] shows that tpi/tp2 
obtained by using the structural data at room temperature is around 1.22 except for 
M =As, in which tpi/tp2 is 1.09. This means that the dimerization gap between the 
SOMO and HOMO bands is small in As, hence resulting in a larger W^somo despite 
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Table 1. List of the transfer energies in the unit of eV for /3-(BDA-TTP)2MF6, where 
LT means a structural data measured in low temperature. 





SbFe (LT) 


PFe 


AsFe 


SbFe 


TaFe 


tpi (eV) 


-0.171 


-0.176 


-0.162 


-0.153 


-0.166 




-0.155 


-0.141 


-0.149 


-0.126 


-0.138 


tql 


-0.070 


-0.065 


-0.066 


-0.055 


-0.061 


tq2 


-0.082 


-0.082 


-0.071 


-0.071 


-0.083 


tc 


0.009 


0.003 


0.008 


0.007 


0.004 


^2c 


0.008 


0.006 


0.005 


0.005 


0.006 


^2p 


0.017 


0.015 


0.018 


0.021 


0.015 


ta 


-0.0002 


-0.002 


0.002 


0.003 


-0.001 


t2q 


0.007 


0.006 


0.007 


0.005 


0.006 


tq3 


0.005 


-0.002 


0.006 


0.003 


-0.002 


tqA 


0.005 


0.006 


0.003 


0.006 


0.009 


trl 


0.014 


0.015 


0.012 


0.014 


0.016 


tr2 


0.015 


0.008 


0.011 


0.008 


0.008 


tpl/tp2 

tql+tq2 


1.10 
0.466 


1.25 
0.464 


1.09 
0.441 


1.21 
0.452 


1.20 
0.474 



the smaller Wsr- This is consistent with the qualitative discussion given in subsection 
13.11 From the viewpoint of the transfer energies, we consider that the difference of 
the anisotropy of the Fermi surface is due to the magnitude of the two- dimensionality, 
{tqi + tq2)/{tpi + tp2)- Our rcsult shows that the average value of {tqi + tq2)/(^pi + ip2) 
over the salts listed in Table[T]is around 0.46, which is smaller than that obtained by the 
extended Hiickel calculation [Sj HI UHl HT] . Therefore, the present Fermi surface exhibits 
an anisotropy in which the a + c-direction is more conductive than in the extended 
Hiickel result. 
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3.3. Random phase approximation and the superconducting gap equation 

We now move on to the RPA analysis for the effective models. To focus on the anion 
dependence, here we concentrate on the models obtained using the lattice parameters at 
room temperature. We show in figur4Zl(a) the spin susceptibility in the Hubbard model 
of /3-(BDA-TTP)2SbF6, where U = 0.38[eV] and the temperature T = 0.008[eV] here. 
The wave vector at which the spin susceptibility is maximized is fc = (0.457r, O.SSvr). 
The arrows in figur^b) corresponds to the peak position of the spin susceptibihty 
in (a), which can be considered as the nesting vector of the Fermi surface. In fact, 
there is a certain overlap between the original Fermi surface and the ones translated by 
these arrows, as shown in figure [7](d). The superconductivity is analysed by solving the 




(C) 
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Figure 7. (a) The calculated spin susceptibility, (b) Arrows correspond to the peak 
position wave vector of the spin susceptibility (nesting vector of the Fermi surface), 
(c) Spin singlet superconducting gap function of the model for /3-(BDA-TTP)2SbF6, 
where U — 0.38, T = 0.008 in 128x128. For this parameter set, the obtained A is 
0.923. (d) The relation between the nesting vector of the Fermi surface and the spin 
susceptibility. 



linearized gap equation that takes into account the pairing interaction mediated by the 
spin fluctuations. The superconducting gap function of the spin-singlet pairing in figure 
W^c) shows that the sign of the gap changes four times along the Fermi surface like in 
a d-wave gap. Considering that the F-B direction corresponds to the a — c direction in 
the lattice structure, the obtained gap function indicates that the hot spots are present 
in the ka ± kc direction. The position of the nodes is different from those obtained 
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in previous theoretical studies that adopt extended Hiickel models (see the appendix) 
[151 [161 [T7] , and hence is not consistent with the STM study [10] . On the other hand, the 
present result is indeed more consistent with the recent measurement of the anisotropy 
of Hc2 in /3-(BDA-TTP)2SbF6 [H]. Further study is necessary to resolve the consistency 
between existing experiments and theories. 

Tc is obtained for the four anions for various on-site interaction U . figure [8] shows 
as functions of the unit-cell volume. There are several features found in the calculation 
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Figure 8. within RPA method for various on-site interaction as a function of the 
unit-ceh volume for ^-(BDA-TTP)2MF6 (Af=P, As, Sb, Ta) of ambient pressure. 



results: (i) Tc is the highest for Sb, (ii) somewhat decreases as we go from P to As, 
and (iii) goes down as we go from Sb to Ta, but it is still finite, close to that of P. 
These features are in accord with the anion dependence of DOS(£^f) or Wsomo shown in 
Fig[ll^c) or (b), but not VTsh or the unit-cell volume. Point (i) is indeed consistent with 
the experimental finding. Point (ii) can also give an understanding for the experimental 
observation that Tc does not increase substantially as we go from P to As. This may at 
first sight be puzzling from the viewpoint of the increase in the unit-cell volume, from 
which one would expect an increase in the density of states. The key for the present 
understanding is the difference in the strength of the dimerization between P and As 
as mentioned previously. As for point (iii), Tc of Ta being somewhat close to P, is not 
consistent with the experiments, where no superconductivity is observed for Ta. The 
discrepancy may be due to effects that are not taken into account in the present model 
(such as phonons), or in the RPA (more sophisticated electron correlation effects). On 
the other hand, the decrease of the Tc with the increase of the unit-cell volume found 
in the present calculation might partially be related to the loss of superconductivity. 
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4. Conclusion 

In the present paper, we have obtained the ab-initio band structure and the effective 
tight-binding model of /3-(BDA-TTP)2MF6 (M =P, As, Sb, Ta), and studied the 
superconducting gap form as well as the anion dependence of Tc assuming spin 
fluctuation mediated pairing. 

The results of ab-initio band calculation for /3-(BDA-TTP)2SbF6 show that the 
direction of the major axis of the cross section of the Fermi surface lies in the F-B 
direction, in contrast to the extended Hiickel result. The density of states at the Fermi 
level is correlated with the band width of SOMO rather than SOMO+HOMO, and hence 
is not necessarily correlated with the unit-cell volume. 

To construct an effective tight-binding model that accurately reproduces the ab- 
initio band structure, we need both nearest and next nearest neighbor transfer energies. 
From the obtained transfer energies, the dimerization of M =As is smaller than that 
of M =P and Sb, and hence the SOMO band width of M =P is smaller than that of 
M =As despite the larger SOMO-^HOMO band width. 

The RPA calculation shows that the wave vector at which the spin susceptibility 
is maximized lies in the a -\- c direction. The gap function of the superconductivity 
mediated by the spin fluctuations changes sign along the Fermi surface like in a ci-wave 
gap, and the hot spots of the gap (gap maxima) is in the ka±kc direction, again different 
from the extended Hiickel model. Tc of M =Sb is the highest among all the anions in 
qualitative agreement with experiments. 

There are some remaining issues that should be resolved in future studies. As 
mentioned above, the relation between the present Fermi surface, those obtained by 
extended Hiickel calculation [3l lU [151 IE]' those obtained in two experiments 
[5l [7] remains as an open question. As for the band calculation, we consider that 
the difference between the present ab-initio result and the extended Hiickel result is 
intimately related to the empirical parameter in the extended Hiickel calculation that 
controls including/excluding the 3d-orbital of the sulfur atoms of the TTP skeleton 
donors [22l [23l [23]. Actually, the extended Hiickel result in which the 3d-orbital of the 
sulfur atoms is excluded [25] is similar to the present ab-initio result. 

Although all experimental and theoretical studies suggest anisotropic supercon- 
ducting gap, the form of the gap function is different between the present and previous 
theoretical studies based on extend Hiickel band [15], [17]. Given the consistency be- 
tween the recent experiment [7] and the present calculation regarding the direction of 
the anisotropy of the Fermi surface, we believe that the gap maxima should lie in the 
ka ± kc direction as far as the spin- fluctuation-mediated pairing mechanism is concerned. 
If the nodes of the superconducting gap do not lie in this direction in the actual materi- 
als, then the present study would conversely suggest that a different pairing mechanism 
may be at work. In this sense, it is highly important to resolve the experimental con- 
troversy regarding the direction of the nodes in the superconducting gap [10|, [TT] . 
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Appendix A. Spin susceptibility and pairing symmetry in the model based 
on the extended Hiickel band 

In this appendix, we show the calculation result obtained by applying the same RPA 
method to the Hubbard model based on the extended Hiickel band. Although the 
present result is essentially the same as those obtained in previous studies [I5l [161 [H] 
here we show the results in the same manner as in figures [7] to make the comparison 
clearer. Here, we concentrate on the model for /3-(BDA-TTP)2SbF6 [3]. The parameter 
values are U = 0.34[eV] and T = 0.008[eV]. 

The calculated spin susceptibility is shown in figure lAl( a). It can be seen that 
it has a broad maximum around the B point, indicating that there is no well defined 
nesting for this Fermi surface. The spin-singlet superconducting gap function in figure 
lAl( b) shows that the sign of the gap changes four times along the Fermi surface, but the 
direction of the nodes differs from that of the ab-initio band model shown in figure [7](c). 
If we consider a wave vector that connects the two gap maxima positions on the Fermi 
surface with different signs of the gap, the vector falls on the broad maximum area in 
the spin susceptibility as shown in figur dAlf c). indicating that the sign change is indeed 
due to the repulsive pairing interaction mediated by the spin fiuctuations. We note that 
the direction of the nodes obtained for this model is in the ka ± kc direction, which is 
in agreement with the STM measurement [10], but not with the in-plane anisotropy of 
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